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Abstract 


This paper assesses the impact of TRMM Microwave Imager (TMI) derived 
surface rainfall data on the numerical simu la tion of Super Typhoon Paka (1997). A series 
of mesoscale numerical simulations is c on ducted for simulating Super Typhoon Paka 
during its mature stage by using the Penn State/NCAR MM5 model with initial 
conditions derived from the Goddard Earth Observing System (GEOS) global analyses 
with and without assimilation of the TMI surface rainfall data. Simulation results clearly 
demonstrate that the GEOS analysis with TMI rainfall data leads to a improved MM5 
simulation of Typhoon Paka in terms of its intensity and kinematical and precipitation 
structures, since inclusion of the rainfall data into initial and boundary condition 
improves the storm environmental conditions (e.g., moisture and circulation). Further 
improvement can also be obtained by incorporating TMI data with mesoscale bogus 
vortex information into the initial conditio ns . 
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1. Introduction 


It has long been recognized that the accurate rainfall measurement over the 
tropical region is one of the most important factors for tropical weather and climate 
studies. Recent numerical studies have also demonstrated that the inclusion of the rainfall 
data into numerical model would impact on the numerical weather forecasts (e.g., 
Krishnamurti et al. 1993; Zupanski and Mesinger 1993). However, the accurate 
measurement of the spatial and temporal variation of tropical rainfall around the globe 
remains one of the critical unsolved problems of meteorology until the lunch of recent 
Tropical Rainfall Measuring Mission (TRMM). During TRMM’s six-year mission, its 
broad sampling footprint between 35°N and 35°S is providing the first detailed and 
comprehensive dataset on the four-dimensional distribution of rainfall and latent heating 
over vastly under sampled oceanic and tropical continental regions. TRMM offers a 
unique opportunity to improve understanding of tropical meteorology and also offers 
great chance to evaluate the impact of rainfall data on the tropical weather forecast. 

Among a lot of the tropical weather events, typhoon is one of the great challenges 
for numerical weather prediction because the conventional observations are usually 
sparse over the tropical ocean. Therefore, it is urgent to incorporate the satellite data into 
numerical weather model. This study assesses the impact of TRMM Microwave Imager 
(TMI) derived surface rainfall data on the numerical simulation of Super Typhoon Paka. 

According to the early study, Super Typhoon Paka was one event frequently 
sampled by TRMM in 1997 (Rodgers et al. 2000). In this study, a set of mesoscale 
numerical simulations is performed forlPaka during its mature stage with the initial 
conditions generated by two different global analyses from the Goddard Earth Observing 
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System (GEOS): one with and one without TRMM rainfall data. The Penn State 
University / National Center for Atmospheric Research (PSU/NCAR) nonhydrostatic 
regional Mesoscale Model version 5 (MM5) is used to conduct the mesoscale 
simulations. Details of the numerical results will demonstrate the impacts of the TRMM 
data on the simulations of Paka, and identify how these impacts are associated with the 
rainfall data. 

Since bogus vortices are often necessary for improving forecasts of mature 
hurricanes (Kurihara et al. 1993; Zou and Xiao 2000; Pu and Braun 2001), as part of this 
evaluation, additional numerical experiments are performed by introducing bogus 
vortices generated by four-dimensional variational data assimilation (4-D VAR) using the 
MM5 adjoint system. In one of the experiments, the bogus vortex is incorporated into the 
mesoscale initial conditions with the TMI data. 

A brief overview of Super Typhoon Paka is given in section 2, and a description 
of the TMI surface rainfall data and GEOS data assimilation are addressed in Section 3. 
The mesoscale model and experimental design are described in section 4, and simulation 
results are given in section 5. A brief summary and discussion are provided in section 6. 

2. A brief overview of Super Typhoon Paka 

Typhoon Paka formed during the first week of December 1997 and underwent 
three periods of rapid intensification over the following two weeks. Rodgers et al. (2000) 
analyzed the detail structural features of Typhoon Paka during its life time using the 
TRMM and the Defense Meteorological Satellite Program’s special sensor 
microwave/imager (SSM/I) satellite observations. Figure 1 shows Paka’s intensity and 
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SSM/I and TMI derived mean total inner- core (within 1 1 1 km of the center) rain rates for 


the period between 9 and 21 Dec 1997. During the earlier periods, which started early on 
December 10, Paka’s maximum wind speed increased from 23 to 60 ms' 1 over a 48-hr 
period with Paka becoming a mature category 3 typhoon on December 1 1 . Later, Paka 
continued to intensify on the following days. By December 18, during the last rapid 
deepening period, Paka became a super typhoon with a maximum wind speed of about 80 
ms' 1 . In order to test the impact of the T MI r ainfall data on the forecast of Paka during its 
mature stage, 0000UTC 12 December 1997 is taken as the initial time for the numerical 
simulations in this study. 

3. Surface rainfall data and GEOS rain fall assimilations 

3.1 TMI Surface Rainfall data 

The surface rainfall information is retrieved from the TMI microwave radiances. 
The basis of the rainfall retrieval method is the “Bayesian” technique described in 
Kummerow et al. (1996) and Olson et a l (1996, 1999). A detailed description of the 
retrieval method utilized in the present study can be found in Olson et al. (1999). The 
retrieval method is supported by a large database of hydrometer profiles that are 
generated from three-dimensional numerical cloud-resolving model simulations (Tao and 
Simpson 1993; Wang et al. 2001). The three-dimensional model fields are used to 
simulate the microwave radiances at the TMI frequencies/polarizations using an 
Eddington radiative transfer method. Given a set of actual TMI observed radiances, the 
retrieved hydrometeor profile is calculated as a weighted mean of the model profiles: 
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only those model profiles associated with radiances similar to the observed radiances are 
designated significant weights. 

3.2 GEOS analysis 

The 6-h averaged surface rainfall estimates derived from the TMT are assimilated 
into the global analysis using a version of GEOS data assimilation system (see Hou et al. 
2000). The analysis system uses the sec ond generation of the GEOS GCM (GEOS-2 
GCM, version 5.9) and version 1.5 of the optimal interpolation (OI) analysis scheme. The 
prognostic variables are potential temperature, specific humidity, surface pressure, and 
winds in the zonal and meridional directions computed at the resolution of 2° latitude, 
2.5° longitude, and 46 a levels from the surface to 4.0 hPa. A unique feature of the GEOS 
data assimilation system is that it uses the “ in cremental analysis update” (IAU) procedure 
of Bloom et al. (1996) to assimilate the rainfall and other observed data. The GEOS 
assimilation is a time-continuous model in te gration with a gradual insertion of IAU “ 
tendencies” into prognostic variables updated from rainfall and other observations every 
6-h. The system eliminates the general spinup problem in rainfall data assimilation. For 
rainfall data assimilation, Hou et al (2000) used a general procedure that minimizes the 
least-square differences between time averaged TMI observations and rain rates 
generated by a column model averaged over a 6-h analysis window. The column model is 
a column version of the GEOS GCM with full model physics, with advection terms as “ 
model forcing” along with the conventional IAU tendencies from a preliminary 6-h 
assimilation. The control variables areT ^tendency correction” of moisture and 
temperature. 
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4. Numerical model and experimental design 


4.1 MM5 model 

The PSU/NCAR MM5 model is used in this study to conduct numerical 
simulations of Paka. The MM5 is a limite d-area, non-hydrostatic primitive equation 
model with multiple options for various physical parameterization schemes. The model 
employs a terrain-following a vertical coordinate, where cr is defined as 
G=(p- pj(p stc - /? top ), p is pressure, and p sfc and p top are the pressures at the surface 

and model top, respectively. Physics options used for this study include the Betts-Miller 
cumulus parameterization, the Goddard c lou d microphysics scheme (Tao and Simpson, 
1993), the Blackadar high-resolution planetary boundary layer parameterization scheme 
(Zhang and Anthes 1982), and the cloud atm ospheric radiation scheme (Dudhia 1993). 
The land surface temperature is predicted using surface energy budget equations as 
described in Grell et al. (1995). For a more detailed description of MM5, the reader is 
referred to Dudhia (1993) and Grell et al. (1995). 

4.2. Experimental design 

A two-way interactive, four-level nested grid technique, is employed to achieve 
the multi-scale simulation. Figure 2 shows the four model domains used and Table 1 
provides specifications for each domain. The outer domains A and B (135 km and 45 km 
horizontal grid spacing) are fixed and are designed to simulate the synoptic-scale and 
mesoscale environment in which the system evolves. The sizes of the domains are chosen 
sufficiently large to minimize the influence of the lateral boundary conditions. The finer 
domains C and D (15 km and 5 km grid spacing) are used to simulate the detailed 
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hurricane-scale flows. The finest domain D i s started at 24 h into the simulation and is 
frequently moved with storm center (e.g., D| t o D 2 ) during the next 36 h simulation. The 
model vertical structure is comprised of 27 <7 levels with the top of the model set at a 
pressure of 50 hPa. The <7 levels are placed at values of 1.0, 0.99, 0.98, 0.96, 0.93, 0.89, 
and then decrease to 0.01 at an interval 0 . 04 . For numerical simulations, model physics 
options are the same for the calculations in all above model domains except there is no 
cumulus parameterization scheme is inc lude d for the calculations in the finest domain 
(domain D, 5 km grid spacing). 

4.3. Model initialization with GEOS analysis 

For the experiments, the initial conditions for the 135-km and 45-km domains 
(domains A and B) are derived from 12-h GEOS analyses. In order to examine the impact 
of TMI data on the storm forecast, two experiments (GEOSO and GEOSTRMM, see 
Table 2) are conducted with initial con dition s generated by two different sets of large- 
scale analyses from the GEOS: a control GEOS analysis data set that does not include 
TMI rainfall data (GEOSO) and a second an alysis data set that does (GEOSTRMM). For 
initialization of grids A and B, the GEOS analysis (without and with the TMI data) 
fields, including potential temperature, specific humidity, surface pressure, and winds in 
the zonal and meridional with the horiz ontal r esolution of 2° latitude and 2.5° longitude 
are interpolated horizontally to the meoscale model grid points. Following the MM5’s 
pre-processing procedure, these interpolated analyses are refined by adding information 
from standard twice-daily rawinsondes and 3-h surface and buoy reports using a Barnes 
objective analysis technique (Manning and Haagenson, 1992). Final analyses are then 
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interpolated to the model a levels. Then, the 15-km domain (domain C) is initialized by 
interpolation (see Grell et al. 1995) of all prognostic variables from the 45-km domain 
using a monotonic interpolation scheme based upon Smolarkiewicz and Grell (1992). 
Domain D is started at 24 h into the forecast and is initialized by interpolation of all 
variables from the 15-km domain. All figures in this paper present results from 15-km 
and 5-km grids. 

Figure 3 shows the MM5 initial con ditions derived from the GEOS analysis 
without and with assimilation of TMI d ata (Fig.3 a and b) in a subset of the 15-km 
domain for fields of sea-level pressure (SLP), wind vectors and divergence at 850 hPa at 
0000UTC 12 December 1997. At the time, Paka was a mature category 3 typhoon, but 
the figures show only weak low pressure system with a broad wind speed maximum to 
the north and east of the circulation center. With assimilation of the TMI rainfall data 
(Fig.3b), only slight differences are foun d i n the SLP, compared with the case without 
TMI data, with the low near 7°N,167°E being about 0.8 hPa deeper. However, The center 
of the low pressure is shifted southward into the confluence circulation (Fig. 3 b and c), 
and the significant differences appear in the wind field with increases of wind speed 
(about 10 ms' 1 ) in the western side and southern side of the confluence center near the 
SLP minimum and therefore the conflue nce bec omes even stronger (Fig.3c). As shown 
in the Figure 3c, the low level divergences are decreased after insertion of TMI data into 
the initial condition. 

Greater differences are observed in the upper level divergences and humidity 
fields. Figure 4 shows the divergence field at lOOhPa and the averaged relative humidity 
between 500 hPa and lOOhPa for both initial conditions without (Fig. 4a) and with (Fig. 
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4b) the TMI data. The figures indicate that assimilation of TMI data leads to notable 
increases of relative humidity and decreases of the divergences in the upper level 
atmosphere (Fig. 4 c). 

4.4. Bogus Vortex 

Since the initial conditions derived from the GEOS analyses contain a poor 
representation of the hurricane vortex, additional numerical experiments are conducted in 
which a bogus vortex is introduced into the initial conditions. 

It has been shown that a reasonable bogus vortex is helpful for hurricane tropical 
cyclone simulations (Kurihara 1993; Leslie et al. 1995). Zou and Xiao (2000) proposed a 
bogus vortex technique that assimilates bogus vortex information using 4-D VAR. The 
method requires two steps: 1) Specification of a bogus vortex by defining the position, 
radius of maximum surface wind (RMW) and minimum SLP of the initial vortex, and 
prescribing a symmetric SLP distribution over the vortex region; 2) assuming that the 
time tendency of SLP is small in a short time period and then assimilating the specified 
bogus SLP field into the numerical mo del with in a 30-minute assimilation window. They 
demonstrated the capability of this technique by using the MM5 adjoint system (Zou et 
al. 1997). Pu and Braun (2000) recently modified this technique to include bogus wind 
information because they found that the wind information was required to ensure gradient 
wind balance, which helps to produce a better vortex structure and prevent vortex spin- 
down problems. In this study, the vorte x tec hnique of Pu and Braun (2000) is used to 
generate the bogus vortex, with assimil ation of the bogus vortex information is applied 
only to the 45-km domain. Based upon the best available estimates, the parameters 
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defining the bogus vortex are the follows: storm centeral pressure p c =955 hPa at the 
typhoon center (7.6°N, 166. 5°E), an environmental pressure of p n = 1012 hPa, maximum 
surface wind V m =51.4 ms" 1 and radius of maximum wind (RMW) R m =135 km. A large 
RMW is used because the 45-km grid spacing is inadequate for resolving the eye for 
realistic vortex of R m . The bogus information extends out to a radius of 350km for both 
SLP and gradient winds. The initial cond i tions derived from 12-h GEOS analyses with 
TMI rainfall, as described in the previous section, are taken as the first guess fields for 
the 4-D VAR system. Thus, the TMI rainfall information and mesosale bogus vortex 
data are incorporated together in the initial conditions. After assimilation of the bogus 
vortex, the 15-km nested domain is initi alized by interpolation (see Grell et al. 1995) of 
all prognostic variables from the 45-km mesh. 

Figure 5 shows the SLP field and w ind vectors and divergence at 850 hPa after 
the bogus vortex assimilation for the 15-km domain. As expected, a strong vortex has 
been successfully introduced into the mesoscale initial fields. The intensity of the bogus 
vortex is close to the estimated intensity of Typhoon Paka. The vortex structure also 
extends throughout the troposphere. Figure 6 compares vertical cross sections of potential 
temperature and wind speed through the typhoon center before and after bogus data 
assimilation, showing that the temperature and wind structure are much improved with 
inclusion of the bogus vortex. A realistic wind structure and strong warm core feature 
near the typhoon center are well resolved in the initial conditions. 

5. Simulation results 

The numerical simulation of Paka is co nducted for 72-h. Figure 7 compares the 
simulated tracks with the observed track for Paka. Results show that inclusion of the 
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TMI data in the initial conditions slows the speed of the storm movement and improves 
the storm track (Table 3). However, this positive improvement does not change the 
direction of movement of the storm so that the track errors still remain large. The bogus 
vortex (BGS) result in significant improvement of the track forecast for Paka, partly 
because it improves the initial location of the vortex, and also because improves the 
interaction of the vortex with the environ men tal flow. 

Figure 8 shows the temporal variation of maximum winds at the lowest model level 
(Fig.8a) and the minimum SLP (or typhoon central pressure, Fig. 8b). The winds are 
compared with the maximum surface wind estimation (courtesy of archives from Unisys 
Corporation on http://weather.uni.sys.com/huiTicane/index.html) . There is no observed 
SLP estimation available for this comparison. The Fig. 8a show substantial improvement 
of the maximum wind forecast when the m o del is started from the initial analysis with the 
TMI data. Although the differences between the “GEOSO” and “GEOSTRMM” are quite 
small at the initial time, significant differences develop after 6 h and even greater 
differences occur after 36 h as the stor m in experiment “GEOSTRMM” becomes a 
typhoon and intensifies rapidly. The storm in the experiment “ GEOSO” undergoes much 
slower intensification and remains a tropical storm until 60-h. The numerical simulation 
with the bogus vortex provides a much better forecast in terms of both maximum wind 
and minimum SLP, comparing well with the estimated maximum wind information 
(Fig.8a). 

Figure 9 compares with the satellite picture of SSM/I brightness temperature from 85 
GHz channel at 0911 UTC 13 December 1997 (Fig.9a) with the forecasted 1-h 
accumulated precipitation at 33 h (corresponding to 0900 UTC 13 December 1997) from 
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the 5-km grid for experiments “GEOSO” (Fig. 9b), “GEOSTRMM” (Fig. 9c) and “BGS” 
(Fig.9d). At the time, Paka was a mature typhoon with a maximum surface wind of about 
53 ms' 1 . Significant differences among the rainfall forecasts are seen between the three 
experiments. Without assimilation of TMI rainfall data (GEOSO), the simulation 
produces only marginal pressure deepening with scattered rainfall on the eastern side of 
the low-pressure system. With assimilation of the TMI rainfall (GEOSTRMM), the 
simulation produces a much more organized outer precipitation structure and a relatively 
strong tropical cyclone. The rainfall pattern reproduces the major features of the satellite 
observed strong convective pattern, but an organized eyewall structure is generally 
lacking. With both TMI rainfall and the bogus vortex (BGS), the simulation produces the 
precipitation features that are consistent with the major observed features, including a 
well-depended eyewall but fails to produce a weak rain band on the western side of the 
storm. 

Figure 10 shows a similar comparison of the rainfall fields, but for the 57-h forecast. 
At this time, experiment “GEOSTRMM” (Fig. 10c) reproduces a mature typhoon with 
central SLP of 974 hPa and maximum surface winds of 45ms' 1 . The rainfall pattern 
matches that observed by SSM/I (Fig. 10a), showing the heavy precipitation on the 
western side of the typhoon eye. Similar features are also observed in experiment “BGS” 
(Fig. lOd) except that the intensity of the storm is even stronger. Without the TMI data 
included in the initial conditions, experiment “GEOSO” remains a tropical storm and the 
precipitation field only shows relatively organized rain bands on the eastern side of the 
storm (Fig. 10b). 
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Figure 1 1 compares the vertical structure of temperature and tangential wind speed 
Field within zonally oriented cross sections through the storm center at 57-h. Similar 
structural features are found in experiments “GEOSTRMM” and “BGS”, as they both 
resolve strong tangential circulation and the warm core feature in the eye very well. On 
the other hand, experiment GEOSO shows a significantly weaker feature. Although the 
wind speeds show a weak maximum ne ar the cen ter of the storm, indicating a weak eye 
wall structure, the overall intensity of the wind speed is much too weak. Only a very 
weak warm temperature anomaly appears in the lower troposphere. 

Since the simulation results with bogus vortex show significant improvement in 
the Paka’s forecast, as a further comparison, an additional experiment (BGSO) is also 
conducted in which a bogus vortex, with the same vortex parameters as in experiment 
BGS, was introduced into the GEOS analysis without assimilation of TMI data. 
Numerical simulation results with experiment BGSO show very similar forecast impacts 
as that in experiment BGS, with storm track and intensity forecasts are improved 
significantly (Fig. 7 and 8). Only slight differences are observed between the both 
experiments. These results might suggest that the bogus vortex can play a dominant role 
in the forecasts of the mature typhoons in terms of their intensity and track. However, 
since the rainfall data are originally assimilated into a global model with coarse 
resolution in this study, it should be necessary to confirm this conclusion in a future 
study, provided the rainfall data are directly assimilated into the mesoscale model. On the 
other hand, as we noted in the figures, the track forecasts in the experiment BGSO are 
slight worse than the forecasts in experiment BGS most of the times. Furthermore, 
distinguishable differences are also observed in the precipitation structure between two 
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experiments during the first 36 h forecast. For instance. Figure 12 shows the 1 h 
accumulated precipitation at 33 h from experiment BGSO, corresponding to the results in 
Figure 9. Comparing two sets of simulations (Fig. 12 and Fig. 9d) with SSM/I 
observations (Fig. 9a), we found that the experiment BGS produces a more realistic 
precipitation structure, which is more consistent with the SSM/I observations. The better 
precipitation structure in BGS is certainly attributed to the impact of including the TMI 
rainfall data in the initial and boundary conditions. The results might also suggest that, 
although the initial vortex can play a dominant role in the forecasts of a matured typhoon, 
assimilation of rainfall data still brings some additional benefits to the forecast, especially 
to the forecast of short-term precipitation. 


6. Summary and Discussion 

Several numerical simulations have been conducted in order to examine the 
impact of assimilated rainfall data on forecast of tropical cyclone structure and intensity. 
Results show significant differences among the first three experiments with initial 
conditions derived from: 1) GEOS analysis without assimilation of TMI rainfall, 2) 
GEOS analysis with assimilation of TMI rainfall data and 3) Similar to 1) and 2) but also 
including mesoscale bogus vortices assimilated into the analysis using 4-D VAR 

technique. The numerical results show that: 

• Assimilation of TMI rainfall dat a into GEOS global analysis results the stronger 
lower level confluences and deeper initial SLP of the storm. It also causes the 
increases of the moisture and divergences in upper level atmosphere and the 
decreases of the divergences in the lower level atmosphere, and therefore lead to 
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the improvement of the environm en tal condition that storm evolutes. As the 
consequent, the storm forecast is improved significantly with the case 
GEOSTRMM in terms of typhoon s tructure and intensity. Since the GEOS large- 
scale analysis does not contain any mesoscale vortex information, experiment 
GEOSTRMM reproduces a storm of typhoon intensity after 36 h forecast, while 
experiment GEOSO requires 60 h of forecast to generate a typhoon. 

• Further forecast improvements are obtained by combining the TMI rainfall data 
with other information, such as a bogus vortex in the initial condition. As shown 
in experiment “BGS”, a bogus vortex that replicates the real intensity of the storm 
can be introduced into the initial conditions. The numerical simulation reproduces 
a very reasonable typhoon track, intensity, and precipitation structure. 

• The bogus vortex can play a dom in a n t role in forecasts of the intensity of mature 
typhoons. However, assimilation of rainfall data into the model still contributes 
additional beneficial impacts, for instance, it helps to produce a more realistic 
short-term precipitation forecast. 

More detail evaluation of TRMM data impact and directly assimilating the TMI 
rainfall into the mesoscale model by 4-D VAR is in progress. Results will be reported in 
the near future. 
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Figure Caption: 

Figure 1. Paka’s intensity (dashed line, ms' 1 ) and SSM/I (open circle) and TMI-derived 

(x) mean total (i.e., combined stratified and convectively generated rain rates) 
inner-core (within 111 km of the center) rain rates (solid line, mm h" 1 ) for the 
period between 9 and 21 Dec 1997. 

Figure 2. Location of the model domains for simulation of Super Typhoon Paka (1997). 
Domain A is the 135-km grid and domain B is the nested, 45-km grid and domain C 
is the 15-km grid used in the forecast. Domain is the 5-km grid and it is started from 
24-h forecast. Domain D moves with the storm center (the domain location may vary 
in different experiment). For experiment “ GEOSTRMM”, domain D move from D| 
to D 2 during 24 to 60-hour forecast. 

Figure 3a. Distributions of the initial sea-level pressure (contour interval, 2 hPa), 
horizontal wind vectors and divergence field (the shaded contour only show the 
divergence less than -0.000007 s' 1 ) at 850 hPa in a part of 15-km domain, for the 
initial conditions in experiment GC OSO (a), experiment GEOSTRMM (b), and the 
differences of all above fields between GEOSTRMM and GEOSO (c). 

Figure 4. The initial divergence (contour interval 2 x 10' 5 s' 1 ) at 150 hPa and averaged 
relative humidity (shaded contour, %) between level 500 hPa to 100 hPa for 
a) GEOSO, b) GEOSTRMM and c) the differences between GEOSTRMM and 
GEOSO. The sign “ X ” denotes the position of the interested circulation center near 
the surface. 

Figure 5. Same as Figure 3a, but for initial condition in experiment “BGS” with the 
contour interval for SLP 8 hPa. 

Figure 6. East-west cross sections of initial potential temperature field (4 K interval) and 
wind speed (5 ms' 1 interval) through the center of storms (corresponding to the 
centers in Fig. 3b and Fig.5) for experiments a) GEOSTRMM and b) BGS. 

Figure 7. Forecasts of the track for all experiments compared to the observed track. 

Center locations along the tracks are indicated every 6 hours. 

Figure 8. Time series (at 6-hour intervals) of a) maximum winds (m s" 1 ) at the lowest 
model level (<j =0.995, approximately 50 m) and b) minimum sea-level pressure 
(hPa). The results during 24h-60h are the forecasts at 5-km grid spacing. 

Figure 9. Comparison of the satellite picture of SSM/I brightness temperature at 85GHz 
Channel at 091 1 UTC 13 December 1997 (a) with the forecasted 1 h accumulated 
precipitation accumulation (contour started from 1 mm h' 1 with interval of 
10 mm h' 1 ) at 33-h forecast (corresponding to 0900 UTC 13 December 1997) at 5- 
km grid spacing, for experiments “GEOSO” (b), “GEOSTRMM” (c) and “BGS” 

(d). 
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Figure 10. Same as Figure 10, except for SSM/I brightness temperature at 0835 UTC 14 
December 1997 (a) and 1 h precipitation accumulation at 57-h forecasts 
(corresponding to 0900 UTC 14 December 1997). 

Figure 1 1 . East-west cross sections of the potential temperature field (4K interval) and 
wind speed (5 m/s interval) through the center of the storms (corresponding to the 
storm centers in Fig 10) at 57-h forecast for experiment a) GEOSO, b) GEOSTRMM 
and c) BGS. 

Figure 12. Same as Fig. 9d, except the forecast from experiment BGS0. 
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Table Caption: 


Table 1. The model design. 


Domain 

Mesh A 

Mesh B 

Mesh C 

Mesh D 

Dimensions 

72x56 

118x88 

202 xl69 

202x169 

Area coverage (km 2 ) 

9720 x 7560 

5310 x 3960 

3030 x 2535 

1010x845 

Grid size (km) 

135 

45 

15 

5 

Time step(s) 

360 

120 

40 

13.3 

Integration hours 

00-72 

00-72 

00-72 

24-60 


Table 2. Experimental design 

Numerical Model initial condition Bogus 

Experiments Vortex 

GEOSO GEOS anlysis without TMI rainfall No 

GEOSTRMM GEOS analysis with TMI rainfall No 

BGS GEOS analysis with TMI and bogus vortex Yes 


BGSO GEOS analysis without TMI rainfall but with bogus vortex Yes 


Table 3. Time series of the 6-h forecasted track error (km) 


Forecast Time (h) 

00 

06 

12 

18 

24 

30 

36 

42 

48 

54 

60 

66 

72 

Track \ 
Error (km)\ 














NumericaT\. \ 
Experiments 

k 













GEOSO 

53 

288 

307 

231 

247 

284 

280 

309 

264 

320 

336 

352 

325 

GEOSTRMM 

41 

156 

210 

218 

211 

211 

227 

215 

160 

185 

174 

192 

170 

BGS 

0 

14 

57 

73 

67 

45 

51 

30 

81 

98 

85 

60 

55 

BGSO 

0 

45 

50 

67 

60 

90 

70 

80 

95 

112 

95 

115 

120 
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